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We perform analytical reductions of one-loop tensor integrals with 5 and 6 legs to scalar master integrals. They 
are based on the use of recurrence relations connecting integrals in different space-time dimensions. The reductions 
are expressed in a compact form in terms of signed minors, and have been implemented in a mathematica package 
called hexagon. m. We present several numerical examples. 



1. Introduction 

Recent years have seen the emergence of first 
results for loop corrections to massive 2 — > 4 scat- 
tering processes |1I2I3) . One of the challenges 
posed by such processes is the need to compute 
one- loop tensor integrals with up to 6 legs. In this 
contribution, we concentrate on integrals with 5 
and 6 legs. 

In 4-dimensional space-time, a linear relation 
must exist between the loop momentum and the 
external momenta of a pentagon integral. Simi- 
larly, if the external momenta of a hexagon are 
4-dimensional, they must be linearly dependent. 
This enables one to reduce hexagons to pen- 
tagons, and pentagons to boxes J4I5I6| . These 
considerations have been extended to dimension- 
ally regularized one-loop tensor integrals [718] . 
Generally, in the reduction of tensor integrals one 
encounters inverse Gram determinants, which can 
vanish at exceptional phase space points and be 
a source of instabilities in a numerical program. 
In the case of 4-dimensional pentagons, it is pos- 
sible to avoid the leading inverse Gram determi- 
nant [S] . Other schemes for the reduction of mul- 
tileg one-loop tensor integrals have been discussed 
inrefs. |10|ll|12|13j . 

Here, we use the methods of ref. [10]. Our 
present goal is to provide compact analytic for- 
mulas for the complete reduction of tensor pen- 
tagons and hexagons to scalar master integrals, 
which are free of leading inverse Gram determi- 



nants. We describe an implementation of these 
formulas in a mathematica package hexagon. m. 

2. Notations 

We consider one-loop, A^-point tensor integrals 
of rank R in d-dimensional space-time. 
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Following Davydychev [Tl], we decompose these 
tensor integrals into a basis of symmetric tensors 
constructed from metric tensors g and the mo- 
menta q.j 
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are Pochhammer symbols 



and the sum runs over non-negative integers such 
that 2X -\- Hi + . . . + UN ~ R- The next step is 
to use recurrence relations to reduce the scalar 



^ In order to be consistent with the conventions of ref. I10| , 
wo have changed the sign of qj as compared with the orig- 
inal definition of ref. |14l . 



coefficients J^^' appearing in the decomposition 
to a set of master integrals. 

It is useful to introduce a notation for certain 
determinants that occur in the recurrence rela- 
tions and their solutions. First, the determinant 
of an {N + 1) X (N + 1) matrix, known as the 
modified Cayley determinant [4] , 
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Although the masses of the propagators appear in 
the coefficients F^ , the determinant Oat does not 
depend on them and it is actually proportional to 
the Gram determinant of the external momenta 
of the A^-point function in eq. ([T]) . All other de- 
terminants we need are signed minors of ()jv, con- 
structed by deleting m rows and m columns from 
Oat, and multiplying with a sign factor. They 
will be denoted by 



Ji 
fci 



J2 
fc2 
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rows ji • • ■ jrn deleted 
columns ki ■ ■ ■ km deleted 
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where sgnr^i and sgnr;.i are the signs of permu- 
tations that sort the deleted rows ji ■ ■ ■ jm and 
columns fci • ■ • km into ascending order. 

Combining integration by parts identities with 
relations connecting integrals in different space- 
time dimensions |15] . one obtains the following 
basic recurrence relations [lOj : 
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where the operator j='= acts by shifting the index 
i^j by ±1. 

3. Pentagons 

A detailed discussion of the second rank pen- 
tagon is given in Ref. [1^. In this section, we 
will consider a third rank tensor integral with in- 
dices 1^1 = . . . = 1/5 = 1, which we write as I^ ^ . 
We assume here and in the next section, that the 
loop momentum k has been shifted in such a way 
that qjv = 0. Applying eq. ([3]) gives integrals in 
space-time dimensions d -I- 4 and d -I- 6 and with 
increased indices. They are reduced back to the 
generic dimension d == 4 — 2e by the recurrence 
relations in eqs. ([T])-®. This involves a division 
by a Gram determinant Oat at each step. The 
leading Gram determinant, Os, can be avoided if 
one is only interested in contractions of the ten- 
sor integral with 4-dimensional objects [S]. This 
is achieved by using the following decomposition 
of the metric tensor, 
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and dropping terms proportional to g^ . After 
further simplifications we obtain: 



jlivX _ 



ir^^ E q^qUkEr^k+Y^g^'^qt^EookAn) 



•■,j,k=l 



fe=l 



where 



with scalar coefficients defined by 



(12) 



(d-5;]^. + i)()^jW(d+2) 



N 



E 
fe=i 



N 



j^'^Hd), 



(8) 



s=l 

+ E 4"^* + E ^l/r^ 



stu 
2 ■> 



(13) 



with 



q4,s _ 
^ijk - 



1 



3(0)5(5^5 

Os 

Ok 

Os 

Os 



■X 



Os 

is 
Oi 
sk 



Os 

5 vJ^'S 
/Os 

5 U^. 



IS 

sk 



+ {i^k) + {j ^ k) 



(14) 




St 



ijk 



3Q5C) 



six 

ust 



5 VJ«/ 5 

Oi\ /wsi 



isf 



fcy 5 \jst 

+ {i^k) + {j ^ k) 



ust 
ist 



(16) 




The decomposition in eq. pip is equivalent to the 
one found in ref. [9J, where the coefhcients Eij^; 
and Eqqj are expressed in terms of tensor 4-point 
functions. Here, instead, they are completely re- 
duced to a basis of scalar master integrals consist- 
ing of boxes /I, vertices /|*, and 2-point functions 
/|*" obtained by removing lines s, s and t, or s, 
t and u from the pentagon. 



4. Hexagons 

If the external momenta of a hexagon are 4- 
dimensional, their Gram determinant vanishes: 
()g = 0, and a linear relation between the propa- 
gators Dj exists: 
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With this relation, any hexagon integral can triv- 
ially be reduced to pentagons. For example, for 
the scalar hexagon, one obtains the well-known 
result ^: 



6 /ON 



-'5 I 



(19) 



where the scalar pentagon /g on the right hand 
side is obtained by removing line r from the 
hexagon /g. In the same way, tensor hexagons 
of rank R can be reduced to tensor pentagons of 
rank R. However, it was noticed in ref. [10' that 
a reduction directly to tensor pentagons of rank 
i? — 1 is also possible: 
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A more general proof of this property was given 
in ref. [13] . By substituting our reduction for- 
mulas for tensor pentagons into eq. (j20p . we can 
immediately express tensor hexagons in terms of 
scalar master integrals. 



5. Numerical results 

The method of reducing six and five point 
Feynman integrals presented in this paper was 
implemented in a MATHEMATICA package cahed 
hexagon. m 17 . It allows to reduce six point in- 
tegrals. Additionally, as hexagons are connected 
with pentagons, the package also allows one to 
reduce pentagons. The present implementation 
includes: 

• six point functions up to rank four 

• five point functions up to rank three 

These tensor ranks are sufficient to get results for 
e.g. NNLO contributions to Bhabha scattering. 
The kinematics used in the package|j is presented 
in Fig. m 

Before using hexagon. m, the package must be 
loaded in a MATHEMATICA environment by execut- 
ing: 

<<hexagon.in 

The package is able to output the full result for a 
six or five point tensor integral, a specific coeffi- 
cient, or a list of all coefhcients for a given rank. 
For a more detailed description, see Table [1] 
hexagon . m is able to generate both analytic and 
numerical results, depending on the user's input. 
It provides coefficients of Lorentz-covariant ten- 
sors, and works in a basis of g'^'^ and internal 
momenta (chords) qi, 



go=0, qn = ^Pt, 



(22) 



see also Fig. [T] In terms of these coefficients, 
the tensor decomposition of pentagons E and 
hexagons F reads: 
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^ The masses and momenta are numbered according to the 
conventions of LoopTools |18| . which are slightly different 
from those used in the previous sections and in ref. 1101 . 




Figure 1 . Momentum flow used in hexagon . m for 
six and five point diagrams. 
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Table 1 

Functions used in the package. 



Six point functions 


Five point 


functions 


RedFO 


scalar 6pt integral 


RedEO 


scalar 5pt integral 


RedFl 


vector 6pt integral 


RedEl 


vector 5pt integral 


RedF2 


rank two 6pt tensor integral 


RedE2 


rank two 5pt tensor integral 


RedF3 


rank three 6pt tensor integral 


RedE3 


rank three 5pt tensor integral 


RedF4 


rank four 6pt tensor integral 






RedFcoef 


coefBcient of given 6pt 


RedEcoef 


coefficient of given 5pt 


RedFget 


all coefficients of given 6pt 


RedEget 


all coefficients of given 5pt 


The basic 


functions have the following arguments, here 


Sij = {Pi+PjY, S,jk = {P^+Pj+Pk)'''■ 


RedFO [pf, 


■ •■ ,pi,Sl2,S23,S34,S45,S56,Si6,Si23,S234,S345 


,m|,...,TO^] 


RedEO[pf, 


• ■ • ,pI, Sl2, S23, S34:, S45, Si5,ml, . 


.,^2] 





For the purpose of checking the correctness of the 
reduction procedures implemented in the pack- 
age, we have made both internal and external 
checks. Internal checks were used for tensor inte- 
grals, and consisted mainly in writing a scalar 
product of internal and external momenta in 
terms of lower rank tensor integrals, which had 
been checked before. External checks were made 
with use of: LoopTools |18j (five point integrals); 
AMBRE [19, with MB.m [20] (five point integrals); 
Sector Decomposition [21j (five and six point 
scalar integrals). 

We present some of the cross-checks in Tabled 
In all these checks, we use LoopTools to calcu- 
late the finite parts of the scalar four, three and 
two point functions which appear after the reduc- 
tion by hexagon. m. We note that, in general, the 
functions defined directly in LoopTools may not 
be sufficient to cover the whole kinematic phase 
space obtained from reduction of six point func- 
tions. In such cases our reduction package must 
be supplemented with additional libraries. 

In the first example in Table [H we use AMBRE 
and MB.m to check the decomposition in the Eu- 
clidean kinematic region, including tensor struc- 
tures. We show the result from hexagon. m for a 
contracted tensor of rank three. At the level typ- 
ical for Monte Carlo calculations, it is in agree- 
ment with the result from Mellin-Barnes integra- 
tion: 0.218885. 

The second example comes from ref. 22], (Ta- 
ble 2, region I). Here, we use the package pi] 
and show for the case of the scalar six point func- 



tion agreement with the calculation using the sec- 
tor decomposition method (typically five digits 
accuracy). In the third example, we have ex- 
tended the scalar case given in [22] (Table 1, re- 
gion III) . Also the tensorial results agree directly 
with LoopTools. 

Finally, for six point tensor integrals, we have 
performed comparisons with a Fortran implemen- 
tation of our reduction formulas by two of us (TD 
and BT), and also with an independent code by 
P. Uwer [23]. Some sample results, for the ran- 
domly chosen phase space point given in table [3l 
are shown in table HI 

6. Summary 

We have described an analytical reduction of 
one-loop tensor integrals with 5 or 6 legs down 
to scalar master integrals, and shown explicitly 
the result for the tensor pentagon of rank three. 
Formulas for other cases and detailed derivations 
will be presented elsewhere. The reduction for- 
mulas have been implemented in a mathematica 
package hexagon . m, which will be made publicly 
available. 
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Table 2 

Numerical cross-checks. 



1) Comparison with AMBRE & MB.m : p'^^PiP^E^.i, 



Point: 

Pi^pI^pI^pI^ 1, pI = 0, mf = to| = 0, m| = to| = m| = 1, 

Sl2 = -3, 523 = -6, S34 = -5, S45 = -7, Sl5 = -2 

In: RedE3 [ pf , . . . ,p^, S12, S23, S34, 345,515,771^, . . . , m^ ] / . {D4->D0,C3->C0,B2->B0} 
Out: 0.218741 



2) Comparison with Sector Decomposition : _Fo 



Point: 

22999? i9 9 9 9 9 9i 

Pi = -P2 = P3 = P4 = -P5 = -Pe = -1> Tni=Tn^= m^ = m| = m^ = m^ = 1, 

512 = 523 = 534 = 545 = 556 = ^16 = ^123 = 5234 = -1, 5345 = -5/2 

In: RedFO [ pf, . . . ,pl, 5i2, 523, 534, 545, 555, 5i6, 5123, 5234, S345, w^, . . . , m^ ] / . {D4->D0} 

Out: 0.013526 

3) Comparison with LoopTools : Eq, Ei, E2, E3, E4, -E'34, -E123, i?oo2 

Point: 

p2 = p2 ^ 0, pI^pI^ 49/256, pI = 9/100, ml = ml = ml = 49/256, ml = ml = 81/1600, 

512 = 4, 523 = ~l/5i 534 = 1/5, 545 = 3/10, 5i5 = -1/2 

In: RedEO[pi, . . . ,^5, 5i2, 523, 534, 545, 515, m^, ...,7715 ]/.D4->D0 

Out: 41.3403 - 45.9721*1 

In: RedEget [rankl ,p\, . . • ,^5, 5i2, 523, 534, 545, 515, mf , . . . , m^ ]/.D4->D0 

Out: eel =-2.38605 + 5.27599*1, ee2 =-5.80875 + 0.597891*1, 

ee3 =-14.4931 + 20.8149*1, ee4 =-11.3362 + 18.1593*1 
In: RedEcoef [ee34 ,Pi, . . . ,p|, 512, 533, 534, 545, 515, to?, . . . , tti^ ] / . {D4->D0,C3->C0} 
Out: 7.1964 + 3.10115*1 

In: RedEcoef [eel23 ,p\,.. • ,pi, S12, S23, S34, 545, 515, to?, . . . , m^ ] / . {D4->D0,C3->C0,B2->B0} 
Out: -0.149527 - 0.31059*1 

In: RedEcoef [ee002 ,p?,.. • ,pi, S12, S23, 534, 545, 515, to?, . . . , m| ] / . {D4->D0,C3->C0,B2->B0} 
Out: 0.154517 - 0.387727*1 
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